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This paper proposes to broaden the canonical formulation of quantum mechanics. Ordinarily, one 
imposes the condition H * — H on the Hamiltonian, where t represents the mathematical operation of 
complex conjugation and matrix transposition. This conventional Hermiticity condition is sufficient 
to ensure that the Hamiltonian H has a real spectrum. However, replacing this mathematical 
condition by the weaker and more physical requirement H* = H, where £ represents combined 
parity reflection and time reversal VT , one obtains new classes of complex Hamiltonians whose 
spectra are still real and positive. This generalization of Hermiticity is investigated using a complex 
deformation H = p 2 + x 2 (ix) e of the harmonic oscillator Hamiltonian, where e is a real parameter. 
The system exhibits two phases: When e > 0, the energy spectrum of H is real and positive as a 
consequence of VT symmetry. However, when — 1 < e < 0, the spectrum contains an infinite number 
of complex eigenvalues and a finite number of real, positive eigenvalues because VT symmetry is 
spontaneously broken. The phase transition that occurs at e — manifests itself in both the 
quantum-mechanical system and the underlying classical system. Similar qualitative features are 
exhibited by complex deformations of other standard real Hamiltonians H — p 2 + x 2N (ix) e with N 
integer and e > — N; each of these complex Hamiltonians exhibits a phase transition at e — 0. These 
"PT-symmetric theories may be viewed as analytic continuations of conventional theories from real 
to complex phase space. 



PACS number(s): 03.65-w, 03.65. Ge, 11.30.Er, 02.60.Lj 
I. INTRODUCTION 

In a recent letter jjl a class of complex quantum- 
mechanical Hamiltonians of the form 



H = p 2 + x 2 (ixf (e real) 



(1.1) 



was investigated. Despite the lack of conventional Her- 
miticity the spectrum of H is real and positive for all 
e > 0. As shown in Fig. |ll| and Fig. 1 of Ref. §, the 
spectrum is discrete and each of the energy levels in- 
creases as a function of increasing e. We will argue below 
that the reality of the spectrum is a consequence of VT 
invariance. 

The operator V represents parity reflection and the 
operator T represents time reversal. These operators are 
defined by their action on the position and momentum 
operators x and p: 



V : x 
T:x 



-x, p 
x, p-> 



-P, 
-p, i -> -i. 



(1.2) 



When the operators x and p are real, the canonical com- 
mutation relation [x,p] = i is invariant under both parity 
reflection and time reversal. We emphasize that this com- 
mutation relation remains invariant under V and T even 
if x and p are complex provided that the above transfor- 
mations hold. In terms of the real and imaginary parts 
of x and p, x — Re x + i Im x and p = Rep + i Imp, we 
have 



V : Hex — > —Rex, Imi — > -Ima;, 
Rep — > —Rep, Imp — > —Imp, 

T : Hex — > Reir, Imi — > -Imo;, 
Rep — > —Rep, Imp —> Imp. 



(1.3) 



While there is as yet no proof that the spectrum of H 
in Eq. (1.1) is real gj, we can gain some insight regarding 
the reality of the spectrum of a 'PT-invariant Hamilto- 
nian H as follows: Note that eigenvalues of the operator 
VT have the form e t6 . To see this, let ^ be an eigenf unc- 
tion of VT with eigenvalue A: VT^ = A 1 ?. Recalling 
that (VT) 2 — 1, we multiply this eigenvalue equation by 
VT and obtain A* A = 1, where we have used the fact that 
i — > —i under VT. Thus, A = e %e . We know that if two 
linear operators commute, they can be simultaneously di- 
agonalized. By assumption, the operator VT commutes 
with H . Of course, the situation here is complicated by 
the nonlinearity of the VT operator (T involves complex 
conjugation). However, let us suppose for now that the 
eigenfunctions ip of H are simultaneously eigenfunctions 
of the operator VT with eigenvalue e 10 . Then applying 
VT to the eigenvalue equation Hip = Eip, we find that 
the energy E is real: E = E* . 

We have numerically verified the supposition that the 
eigenfunctions of H in Eq. (1.1) arc also eigenfunctions 
of the operator VT when e > 0. However, when e < 0, 
the VT symmetry of the Hamiltonian is spontaneously 
broken; even though VT commutes with H , the eigen- 
functions of H are not all simultaneously eigenfunctions 
of VT . For these eigenfunctions of H the energies are 
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complex. Thus, a transition occurs at e = 0. As e goes 
below 0, the eigenvalues as functions of e pair off and be- 
come complex, starting with the highest-energy eigenval- 
ues. As e decreases, there are fewer and fewer real eigen- 
values and below approximately e = —0.57793 only one 
real energy remains. This energy then begins to increase 
with decreasing e and becomes infinite as e approaches 
— 1. In summary, the theory defined by Eq. ( |l.l| ) exhibits 
two phases, an unbroken-symmetry phase with a purely 
real energy spectrum when e > and a spontaneously- 
broken-symmetry phase with a partly real and partly 
complex spectrum when e < 0. 

A primary objective of this paper is to analyze the 
phase transition at e = 0. We will demonstrate that this 
transition occurs in the classical as well as in the quan- 
tum theory. As a classical theory, the Hamiltonian H de- 
scribes a particle subject to complex forces, and therefore 
the trajectory of the particle lies in the complex- a; plane. 
The position and momentum coordinates of the particle 
are complex functions of t, a real time parameter. We 
are interested only in solutions to the classical equations 
of motion for which the energy of the particle is real. We 
will see that in the PT-symmetric phase of the theory, 
the classical motion is periodic and is thus a complex gen- 
eralization of a pendulum. We actually observe two kinds 
of closed classical orbits, one in which the particle oscil- 
lates between two complex turning points and another 
in which the particle follows a closed orbit. In many 
cases these closed orbits lie on an elaborate multishected 
Riemann surface. On such Riemann surfaces the closed 
periodic orbits exhibit remarkable knot-like topological 
structures. All of these orbits exhibit PT symmetry; 
they are left-right symmetric with respect to reflections 



about the imaginary- a; axis in accordance with Eq. (1.3) 



In the broken-symmetry phase classical trajectories are 
no longer closed. Instead, the classical path spirals out 
to infinity. These spirals lack PT symmetry. 

There have been many previous instances of non- 
Hermitian PT-invariant Hamiltonians in physics. Ener- 
gies of solitons on a complex Toda lattice have been found 
to be real fl. Hamiltonians rendered non-Hermitian by 
an imaginary external field have been used to study popu- 
lation biology j|| and to study derealization transitions 
such as vortex flux-line depinning in type-II supercon- 
ductors H. In these cases, initially real eigenvalues bi- 
furcate into the complex plane due to the increasing ex- 
ternal field, indicating the growth of populations or the 
unbinding of vortices. 

The PT-symmetric Hamiltonian considered in this pa- 
per has many generalizations: (i) Introducing a mass 



term of the form 



yields a theory that exhibits sev- 



eral phase transitions; transitions occur at e = — 1 and 
e = — 2 as well as at e = |l|. (ii) Replacing the condition 
of Hermiticity by the weaker constraint of PT-symmetry 
also allows one to construct new classes of quasi-exactly 
solvable quantum theories B . (iii) In this paper we con- 



sider complex deformations of real Hamiltonians other 
than the harmonic oscillator. We show that Hamiltoni- 
ans of the form 



H =p 2 + x 2K 



(i X y 



(1.4) 



have the same qualitative properties as H in Eq. (1.1). As 
e decreases below 0, all of these theories exhibit a phase 
transition from an unbroken PT-symmetric regime to a 
regime in which PT symmetry is spontaneously broken. 

The Hamiltonian H in ([hi]) is especially interesting 
because it can be generalized to quantum field theory. A 
number of such generalizations have recently been exam- 
ined. The PT-symmetric scalar field theory described by 
the Lagrangian |7j 

£= i(90) 2 + imV+50 2 W) e (e>0) (1.5) 

is intriguing because it is not invariant under parity re- 
flection. This is manifested by a nonzero value of {<j)). 
It is interesting that this broken symmetry persists even 
when e > is an even integer Q. The Hamiltonian for 
this theory is not Hermitian and, therefore, the theory is 
not unitary in the conventional sense. However, there is 
strong evidence that the spectrum for this theory is real 
and bounded below. For e = 1 one can understand the 
positivity of the spectrum in terms of summability. The 
weak-coupling expansion for a conventional g^ 3, theory 
is real, and apart from a possible overall factor of g, the 
Green's functions are formal power series in g 1 . These 
series are not Borel summable because they do not alter- 
nate in sign. Nonsummability reflects the fact that the 
spectrum of the underlying theory is not bounded below. 
However, when we replace g by ig, the perturbation se- 
ries remains real but now alternates in sign. Thus, the 
perturbation series becomes summable, and this suggests 
that the underlying theory has a real positive spectrum. 

Replacing conventional gef) 4 or gcf> 3 theories by PT- 
symmetric —gcf> 4 or igtf) 3 theories has the effect of re- 
versing signs in the beta function. Thus, theories that 
are not asymptotically free become asymptotically free 
and theories that lack stable critical points develop such 
points. There is evidence that —gcf> 4 in four dimensions 
is nontrivial ||. 

Supersymmetric quantum field theory that is PT in- 
variant has also been studied g). When we construct 
a two-dimensional supersymmetric quantum field theory 
by using a superpotential of the form S(4>) = —ig{i4>) 1+e , 
the supersymmetric Lagrangian resulting from this su- 
perpotential is 

£ = \ W) 2 + \$W> + \s'(4>)^ + \W)] 2 
= \ W) 2 + + \g(i + e)(i<t>) e W> 
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\2+2e 



(1.6) 
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where ip is a Majorana spinor. The Lagrangian (1.3) has 
a broken parity symmetry. This poses the question, Does 
the parity violation induce a breaking of supersymmetry? 
To answer this question, both the ground-state energy 
Eq and the fermion-boson mass ratio R were calculated 
as series in powers of the parameter e. Through second 
order in e, Eq = and R = 1, which strongly suggests 
that supersymmetry remains unbroken. We believe that 
these results are valid to all orders in powers of e. This 
work and our unpublished numerical studies of SUSY 
quantum mechanics show that complex deformations do 
not break supersymmetry. 

Quantum field theories having the property of VT in- 
variance exhibit other interesting features. For example, 
the Ising limit of a 'PT-invariant scalar quantum field 
theory is intriguing because it is dominated by solitons 
rather than by instantons as in a conventional quantum 
field theory |Tc|| . In addition, a model of 'PT-invariant 
quantum electrodynamics has been studied [p"l| . The 
massless theory exhibits a stable, nontrivial fixed point at 
which the renormalized theory is finite. Moreover, such a 
theory allows one to revive successfully the original elec- 
tron model of Casimir. 

Since (4>) =/= in 'PT-symmetric theories, one can in 
principle calculate directly (using the Schwinger-Dyson 
equations, for example) the real positive Higgs mass in 
a renormalizable PT-symmetric theory in which symme- 
try breaking occurs naturally. No symmetry-breaking 
parameter needs to be introduced. This most intriguing 
idea could lead to an experimental vindication of our pro- 
posed generalization of the notion of Hermiticity to PT 
symmetry. 

This paper is organized as follows: In Sec. p| we study 
the classical version of the Hamiltonian in Eq. fll.ip . The 
behavior of classical orbits reveals the nature of the phase 
transition at e = 0. Next, in Sec. Ill we analyze the 
quantum version of this Hamiltonian. We derive several 
asymptotic results regarding the behavior of the energy 
levels near the phase transition. In Sec. we discuss 
the classical and quantum properties o f the broad class 
of PT-symmetric Hamiltonians in Eq. (|l.4| ) of which H 
in Eq. (LI) is a special case. Finally, in Sec. we study 
complex deformations of nonanalytic potentials. 



II. CLASSICAL THEORY 

The classical equa tion of motion for a particle de- 
scribed by H in (LI) is obtained from Hamilton's equa- 
tions: 



dx dH 

dt dp 
dp dH 
dt dx 

Combining these two equations gives 



2p, 

= i(2 + e){ix) 1+e . 



(2.1) 



fx 
~dP 



= 2i{2 + e){ixf 



(2.2) 



which is the complex version of Newton's second law, 
F = ma. 

Equation (2.2) can be integrated once to give ]L2| 



1 dx 
2~dt 



±^/E + (ix) 2 +c, 



(2.3) 



where E is the energy of the classical particle (the time- 
independent value of H ) . We treat time t as a real vari- 
able that parameterizes the complex path x(t) of this 
particle. 

This section is devoted to studying and classifying the 
solutions to Eq. (2.3). By virtue of the PT invariance 
of the Hamiltonian H, it seems reasonable to restrict 
our attention to real values of E. Given this restriction, 
we can always rescale x and t by real numbers so that 
without loss of generality Eq. (12. 31) reduces to 



dx 
~dt 



±y/l + (ix) 



2+e 



(2.4) 



The trajectories satisfying Eq. (2.4) lie on a multi- 
shected Ricmann surface. On this surface the function 
v/l + (ix) 2+£ is single- valued. There are two sets of 
branch cuts. The cuts in the first set radiate outward 
from the roots of 



1 + ( M ) 2 + e = 



(2.5) 



These roots are the classical turning points of the motion. 
There are many turning points, all lying at a distance of 
unity from the origin. The angular separation between 
consecutive turning points is 2it/(2 + e). The second set 
of branch cuts is present only when e is noninteger. In 
order to maintain explicit PT symmetry (left-right sym- 
metry in the complex- x plane), we choose these branch 
cuts to run from the origin to infinity along the positive 
imaginary axis. 



A. Case e = 



Because the classical solutions to Eq. (2.4) have a very 
elaborate structure, we begin by considering some special 
values of e. The simplest case is e = 0. For this case there 
are only two turning points and these lie on the real axis 
at ±1. 



In order to solve Eq. (2.4) we need to specify an initial 
condition x(0). The simplest choice for x(0) is a turning 
point. If the path begins at ±1, there is a unique direc- 
tion in the complex- x plane along which the phases of 
the left side and the right side of Eq. (|]4|) agree. This 
gives rise to a trajectory on the real axis that oscillates 
between the two turning points. This is the well-known 
sinusoidal motion of the harmonic oscillator. 
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Note that once the turning points have been fixed the 
energy is determined. Thus, choosing the initial position 
of the particle determines the initial velocity (up to a 
plus or minus sign) as well. So, if the path of the parti- 
cle begins anywhere on the real axis between the turning 
points, the initial velocity is fixed up to a sign and the 
trajectory of the particle still oscillates between the turn- 
ing points. 

Ordinarily, in conventional classical mechanics the only 
possible initial positions for the particle lie on the real- a; 
axis between the turning points because the velocity is 
real; all other points on the real axis lie in the classically 
forbidden region. However, because we are analytically 
continuing classical mechanics into the complex plane, 
we can choose any point x(0) in the complex plane as an 
initial position. For all complex initial positions outside 
of the conventional classically allowed region the classical 
trajectory is an ellipse whose foci are the turning points. 
The ellipses are nested because no trajectories may cross. 
(See Fig. 0.) The exact solution to Eq. (0) is 

x(t) — cos[arccosx(0) ± t], (2-6) 

where the sign of t determines the direction (clockwise 
or anticlockwise) in which the particle traces the ellipse. 
For any ellipse the period of the motion is 2ir. The pe- 
riod is the same for all trajectories because we can join 
the square-root branch cuts emanating from the turning 
points, creating a single finite branch cut lying along the 
real axis from x = —1 to x = X. The complex path inte- 
gral that determines the period can then be shrunk (by 
Cauchy's theorem) to the usual real integral joining the 
turning points. 

2.5 | , 1 , , , 1 




-2.5 1 1 1 1 1 1 1 
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FIG. 1. Classical trajectories in the complex- a; plane for the 
harmonic oscillator whose Hamiltonian is H — p 2 + x 2 . These 
trajectories represent the possible paths of a particle whose 
energy is E = 1. The trajectories are nested ellipses with 
foci located at the turning points at x — ±1. The real line 
segment (degenerate ellipse) connecting the turning points is 
the usual periodic classical solution to the harmonic oscillator. 
All closed paths [see Eq. (2.6)] have the same period 2n. 



Finally, we remark that all of the classical paths (el- 
liptical orbits) are symmetric with respect to parity V 
(reflections through the origin) and time reversal T (re- 
flections about the real axis), as well as VT (reflections 
about the imaginary axis). Furthermore, V and T indi- 
vidually preserve the directions in which the ellipses are 
traversed. 

B. Case e = 1 

The case e = 1 is significantly more complicated. Now 
there are three turning points. Two are located below 
the real axis and these are symmetric with respect to the 
imaginary axis: = e~ 5j7r / 6 and x + = e~ 47r / 6 . That is, 
under VT reflection x~ and x + are interchanged. The 
third turning point lies on the imaginary axis at xq = i. 

As in the case e = 0, the trajectory of a particle that 
begins at the turning point X- follows a unique path in 
the complex- x plane to the turning point at x+. Then, 
the particle retraces its path back to the turning point 
at X-, and it continues to oscillate between these two 
turning points. This path is shown on Fig. ||. The period 
of this motion is 2-\/37rr(|)/r(|). The periodic motion 
between x± is clearly time-reversal symmetric. 




-6 1 1 1 1 1 1 1 1 1 

-4-3-2-101234 

FIG. 2. Classical trajectories in the complex- £ plane for 
a particle described by the Hamiltonian H = p 2 + ix 3 and 
having energy E = 1. An oscillatory trajectory connects the 
turning points x±. This trajectory is enclosed by a set of 
closed, nested paths that fill the finite complex- a; plane ex- 
cept for points on the imaginary axis at or above the turning 
point xo — i. Trajectories originating at one of these excep- 
tional points go off to ioo or else they approach xo, stop, turn 
around, and then move up the imaginary axis to ioo. 

A particle beginning at the third turning point xq 
exhibits a completely distinct motion: It travels up 
the imaginary axis and reaches ioo in a finite time 
V / 7rr(|)/r(|). This motion is not periodic and is not 
symmetric under time reversal. 

Paths originating from all other points in the finite 
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complex- a; plane follow closed periodic orbits. No two 
orbits may intersect; rather they are all nested, like the 
ellipses for the case e = 0. All of these orbits encircle the 
turning points x± and, by virtue of Cauchy's theorem, 
have the same period 2v / 37rr(|)/r(|) as the oscillatory 
path connecting x±. Because these orbits must avoid 
crossing the trajectory that runs up the positive imagi- 
nary axis from the turning point xq — i, they are pinched 
in the region just below xo, as shown on Fig. |[ 

As these orbits become larger they develop sharper in- 
dentations in the vicinity of xq. We observe that the 
characteristic radius of a large orbit approaches the re- 
ciprocal of the distance d between xo and the point where 
the orbit intersects the positive imaginary axis. Thus, it 
is appropriate to study these orbits from the point of 
view of the renormalization group: We scale the distance 
d down by a factor L and then plot the resulting orbit on 
a graph whose axis are scaled down by the same factor 
L. Repeated scaling gives a limiting orbit whose shape 
resembles a cardioid (see Fig. ||). The equation of this 
limiting orbit is obtained in the asymptotic regime where 
we neglect the dimensionless energy 1 in Eq. (2.4): 



d it =±[lxf ' 2 - 



(2.7) 



The solution to this differential equation, scaled so that 
it crosses the negative imaginary axis at —3?, is 



x(t) = 



At 



(t + 2i/Vsy 



(— oo < t < oo). 



(2.8) 



This curve is shown as the solid line in Fig. ||. (Strictly 
speaking, this curve is not a true cardioid, but its shape 
so closely resembles a true cardioid that we shall refer to 
it in this paper as the limiting cardioid.) 

In the infinite scaling limit all periodic orbits [all 
these orbits have period 2V37rT(|)/r(|)], which origi- 
nally filled the entire finite complex- a; plane, have been 
squeezed into the region inside the limiting cardioid ( 2 



The nonperiodic orbit still runs up the positive imaginary 
axis. The obvious question is, What complex classical dy- 
namics is associated with all of the other points in the 
scaled complex- a; plane that lie outside of the limiting 
cardioid? We emphasize that all of these points were 
originally at infinity in the unsealed complex- a; plane. 

We do not know the exact answer to this question, 
but we can draw a striking and suggestive analogy with 
some previously published work. It is generally true that 
the region of convergence in the complex- a; plane for an 
infinitely iterated function is a cardioid-shaped region. 
For example, consider the continued exponential function 



f(x) = e" 



(2.9) 



The sequence e x , e xe , • • • is known to converge in 
a cardioid-shaped region of the complex- a; plane (see 



Figs. 2-4 in Ref. |13|). It diverges on the straight line 
that emerges from the indentation of the cardioid. The 
remaining part of the complex- a; plane is divided into an 
extremely elaborate mosaic of regions in which this se- 
quence converges to limit cycles of period 2, period 3, pe- 
riod 4, and so on. These regions have fractal structure. It 
would be interesting if unbounded complex classical mo- 
tion exhibits this remarkable fractal structure. In other 
words, does the breaking of V and T symmetry allow for 
unbounded chaotic solutions? 
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FIG. 3. Approach to the limiting cardioid in Eq. (2.8). As 
the orbits shown in Fig. ^approach the turning point xq, they 
get larger. Using a renormalization-group approach, we plot 
successively larger orbits (one such orbit is shown as a dashed 
line) scaled down by the characteristic size of the orbit. The 
limiting cardioid is indicated by a solid line. The indentation 
in the limiting cardioid develops because classical trajectories 
may not intersect and thus must avoid crossing the trajectory 
(shown in Fig. ^ on the imaginary axis above xq. 



C. Case e = 2 

When e = 2 there are four turning points, two located 
below the real axis and symmetric with respect to the 
imaginary axis, x\ = e~ 3l7T / 4 and X2 — e~ 47r / 4 , and two 
more located above the real axis and symmetric with re- 
spect to the imaginary axis, x% = e 47r / 4 and x^ — e 3l7T / 4 . 
Classical trajectories that oscillate between the pair x\ 
and X2 and the pair X3 and X4 are shown on Fig. ^. The 
period of these oscillations is 2v / 27rr(|)/r(|). Trajecto- 
ries that begin elsewhere in the complex- a; plane are also 
shown on Fig. ||. Note that by virtue of Cauchy's theo- 
rem all these nested nonintersecting trajectories have the 
same period. All motion is periodic except for trajecto- 
ries that begin on the real axis; a particle that begins on 
the real- a; axis runs off to ±00, depending on the sign of 
the initial velocity. These are the only trajectories that 
are nonperiodic. 

The rescaling argument that gives the cardioid for the 
case e = 1 yields a doubly-indented cardioid for the case 
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e = 2 (see Fig 
Fig. 5 of Ref. fl3fl 



. This cardioid is similar to that in 
However, for the case e = 2 the limit- 
ing double cardioid consists of two perfect circles, which 
are tangent to one another at the origin x = 0. Circles 
appear because at e = 2 in the scaling limit the equation 
corresponding to ( |2.7D is ^ = ±x 2 . The solutions to this 
equation are the inversions x(t) — ifrj, which map the 
real-i axis into circles in the complex- x plane. 



One other class of trajectory is possible. If the initial 
position of the classical particle lies on the imaginary 
axis at or above the turning point at z, then depending 
on the sign of the initial velocity, the particle either runs 
off to zoo or it approaches the turning point, reverses its 
direction, and then goes off to zoo. These purely imagi- 
nary paths are the only possible nonperiodic trajectories. 
They are also shown on Fig. ^|. 
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FIG. 4. Classical trajectories in the complex- 1 
a particle described by the Hamiltonian H — p 2 
having energy £7 = 1. There are two oscillatory trajecto- 
ries connecting the pairs of turning points xi and X2 in the 
lower-half x-plane and x$ and x<t in the upper-half s-plane. 
[A trajectory joining any other pair of turning points is for- 
bidden because it would violate VT (left-right) symmetry.] 
The oscillatory trajectories are surrounded by closed orbits 
of the same period. In contrast to these periodic orbits there 
is a class of trajectories having unbounded path length and 
running along the real-x axis. These are the only paths that 
violate time-reversal symmetry. 
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FIG. 5. Limiting double cardioid for the case e = 2. As the 
orbits in Fig. ^| approach the real axis, they get larger. If we 
scale successively larger orbits down by their characteristic 
size, then in the limiting case the orbits approach two circles 
tangent at the origin. In this limit the four turning points in 
Fig. coalesce at the point of tangency. 



D. Case e = 5 

When e = 5 there are seven turning points, one located 
at i and three pairs, each pair symmetric with respect to 
reflection about the imaginary axis (VT symmetric) . We 
find that each of these pairs of turning points is joined by 
an oscillatory classical trajectory. (A trajectory joining 
any other two turning points would violate VT symme- 
try) Surrounding each of the oscillatory trajectories are 
nested closed loops, each loop having the same period 
as the oscillatory trajectory it encloses. These classical 
trajectories are shown on Fig. |^. The periods for these 
three families of trajectories are 

where 9 — bn/lA for the lowest pair of turning points, 
9 = tt/14 for the middle pair, and 9 = 3ir/14 for the pair 
above the real axis. 
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FIG. 6. Classical trajectories in the complex- x plane for 
a particle described by the Hamiltonian H = p 2 + ix 7 and 
having energy £7 = 1. Shown are oscillatory trajectories sur- 
rounded by periodic trajectories. Unbounded trajectories run 
along the positive-imaginary axis above x = i. 
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E. General case: Noninteger values of e > 



Because Eq. ( 2A ) contains a square-root fun ctio n, the 
turning points, which are solutions to Eq. ( |2.5|) , are 
square-root branch points for all values of e. Thus, in 
principle, the complex trajectories x(t) lie on a multi- 
sheeted Riemann surface. However, when e is a nonneg- 
ative integer, we can define the branch cuts so that the 
classical trajectories satisfying Eq. ( |2.4| ) never leave the 
principal sheet of this Riemann surface. We do this as 
follows: We choose to join the 'PT-symmetric (left-right- 
symmetric) pairs of turning points by branch cuts that 
follow exactly the oscillatory solutions connecting these 
pairs. (There are three such pairs in Fig. 0, two in Fig. ^, 
and one in Figs. || and [l].) If e is odd, there is one ex- 
tra turning point that lies on the positive imaginary axis 
(see Figs, g and ^); the branch cut emanating from this 
turning point runs up the imaginary-x axis to ioo. Since 
classical paths never cross, there are no trajectories that 
leave the principal sheet of the Riemann surface. 

When e is noninteger, we can se e fro m the argument 
of the square-root function in Eq. (2^4) that there is an 
entirely new branch cut, which emerges from the origin 
in the complex- a; plane. To preserve VT symmetry we 
choose this branch cut to run off to oo along the positive- 
imaginary x-axis. If e is rational, the Riemann surface 
has a finite number of sheets, but if e is irrational, then 
there are an infinite number of sheets. 

If a classical trajectory crosses the branch cut ema- 
nating from the origin, then this trajectory leaves the 
principal sheet of the Riemann surface. In Fig. [t] we il- 
lustrate some of the possible classical trajectories for the 
case e = tt — 2. The top plot shows some trajectories 
that do not cross the positive-imaginary x-axis and thus 
do not leave the principal sheet of the Riemann surface. 
The trajectories shown are qualitatively similar to those 
in Fig. |2|; all trajectories have the same period. 

In the middle plot of Fig. [?| is a trajectory that crosses 
the positive-imaginary x-axis and visits three sheets of 
the Riemann surface. The solid line and the dotted line 
outside of the solid line lie on the principal sheet, while 
the remaining two portions of the dotted line lie on two 
other sheets. Note that this trajectory does not cross 
itself; we have plotted the projection of the trajectory 
onto the principal sheet. The trajectory continues to 
exhibit VT symmetry. The period of the trajectory is 
greater than that of the period of the trajectories shown 
in the top plot. This is because the trajectory encloses 
turning points that are not on the principal sheet. In 
general, as the size of the trajectory increases, it encloses 
more and more complex turning points; each time a new 
pair of turning points is surrounded by the trajectory the 
period jumps by a discrete quantity. 

Although the trajectory in the bottom plot in Fig. |^ 
has the same topology as that in the middle plot, it is 
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larger. As the trajectory continues to grow, we observe 
a phenomenon that seems to be universal; namely, the 
appearance of a limiting cardioid shape (solid line) on 
the principal surface. The remaining portion of the tra- 
jectory (dotted line) shrinks relative to the cardioid and 
becomes compact and knot-like. 
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FIG. 7. Classical trajectories for H — p 2 — (ix)* corre- 
sponding to the case e = n — 2. Observe that as the classical 
trajectory increases in size, a limiting cardioid appears on the 
principal sheet of the Riemann surface. On the other sheets 
the trajectory becomes relatively small and knot-like. 
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In Fig. |]we examine the case e = 0.5. In this figure we 
observe behavior that is qualitatively similar to that seen 
in Fig. fj]; namely, as the trajectory on the principal sheet 
of the Riemann surface becomes larger and approaches a 
limiting cardioid, the remaining portion of the trajectory 
becomes relatively small and knot-like. 
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FIG. 8. Classical trajectories for the case e = 0.5. As the 
classical path on the principal sheet of the Riemann surface 
increases in size it approaches a limiting cardioid, just as in 
Fig. (?]. The remaining portion of the path becomes relatively 
small and knot-like. 



To summarize, for any e > the classical paths are 
always VT symmetric. The simplest such path describes 
oscillatory motion between the pair of turning points that 
lie just below the real axis on the principal sheet. In 
general, the period of this motion as a function of e is 
given by 



r(gf) 



67T 



2c 



(2.10) 



Other closed paths having more complicated topologies 
^and longer periods) also exist, as shown in Figs. |t] and 

Whenever the classical motion is periodic, we expect 
the quantized version of the theory to exhibit real eigen- 
values. Although we have not yet done so, we intend to 
investigate the consequences of quantizing a theory whose 
underlying classical paths have complicated topological 
structures traversing several sheets of a Riemann surface. 
The properties of such a theory of quantum knots might 
well be novel. 



F. Case -1 < e < 

Classical paths for negative values of e are fundamen- 
tally different from those corresponding to nonnegative 
values of e; such paths no longer exhibit VT symmetry. 
Furthermore, we no longer see paths that are periodic; all 
paths eventually spiral outwards to infinity. In general, 
the time that it takes for a particle to reach infinity is 
infinite. 

We interpret the abrupt change in the global nature of 
the classical behavior that occurs as e passes through 
as a cha nge in phase. For all values of e the Hamiltonian 
in Eq. (1.1) is VT (left-right) symmetric. However, for 
e < the solutions cease to exhibit VT symmetry. Thus, 
we say that e > is a PT-symmetric phase and that 
e < is a spontaneously broken 'PT-symmetric phase. 

To illustrate the loss of VT (left-right) symmetry, we 
plot in Fig. ^| the classical trajectory for a particle that 
starts at a turning point x_ = — 7r in the second 
quadrant of the complex- a; plane (Rex < 0, Imi > 0) 
for three values of e: —0.2, —0.15, and —0.1. We ob- 
serve that a path starting at this turning point moves to- 
ward but misses the PT-symmetric turning point x + = 
—TTjf^- because it crosses the branch cut on the positive- 



4+2e 

imaginary x-axis. This path spirals outward, crossing 
from sheet to sheet on the Riemann surface, and even- 
tually veers off to infinity asymptotic to the angle 6oo, 
where 

#oo = -^7^. (2.11) 

This formula shows that the total angular rotation of 
the spiral is finite for all e 7^ but becomes infinite as 
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e — > _ . In the top figure (e = —0.2) the spiral makes 2 j 
turns before moving off to infinity; in the middle figure 
(e = —0.15) the spiral makes 3^ turns; in the bottom 

| turns. 



figure (e = —0.1) the spiral makes 4 ■ 




-30 



-70 -60 -50 -40 -30 -20 -10 



10 



20 



40 
20 

-20 
-40 
-60 
-80 
-100 
-120 

-30 -20 -10 10 20 30 40 50 60 

FIG. 9. Classical trajectories that violate VT symmetry. 
The top plot corresponds to the case e = —0.2, the middle plot 
to £ = —0.15, and the bottom plot to e = —0.1. The paths 
in each plot begin at a turning point and spiral outwards to 
infinity in an infinite amount of time. 

Note that the spirals in Fig. pass many classical 




turning points as they spiral clockwise from X-. [From 
Eq. ( |2.5| ) we see that the nth turning point lies at the an- 
gle 4 " K (%- corresponds ton — 0).] As e approaches 
from below, when the classical trajectory passes a new 
turning point, there is a corresponding merging of the 
quantum energy levels as shown in Fig. 11. As pointed 
out in Rcf. 0| , this correspondence becomes exact in the 
limit e — * 0~ and is a manifestation of Ehrenfest's theo- 
rem. 



G. Case e 



For this special case we can solve the equation (2.4) 
exactly. The result, 



x(t) = I 1 - Y + -r \ i + bt (b real), (2.12) 

represents a family of parabolas that are symmetric with 
respect to the imaginary axis (see Fig. [l0|). Note that 
there is one degenerate parabola corresponding to b = 
that lies on the positive imaginary axis above i. 




FIG. 10. Classical trajectories in the complex- a; plane for a 
particle described by the Hamiltonian H = p' 2 — ix and having 
energy E = 1. Shown are parabolic trajectories and a turning 
point at i. All trajectories are unbounded. 



III. QUANTUM THEORY 



In this section we discuss the quantum properties of 



the Hamiltonian H in Eq. (1.1). The spectrum of this 
Hamiltonian is obtained by solving the corresponding 
Schrodinger equation 



ip"(x) + [x 2 (ix) e -E\tp(x) = 



(3.1) 



subject to appropriate boundary conditions imposed in 
the complex- a; plane. These boundary conditions are de- 
scribed in Ref. [Q . A plot of the spectrum of H is shown 
in Fig. |ll|. 
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FIG. 11. Energy levels of the 

Hamiltonian H — p 2 + x 2 {ix) e as a function of the param- 
eter e. There are three regions: When e > 0, the spectrum 
is real and positive and the energy levels rise with increasing 
e. The lower bound of this region, e = 0, corresponds to the 
harmonic oscillator, whose energy levels are E n = 2n + 1. 
When — 1 < e < 0, there are a finite number of real positive 
eigenvalues and an infinite number of complex conjugate pairs 
of eigenvalues. As e decreases from to — 1, the number of 
real eigenvalues decreases; when e < —0.57793, the only real 
eigenvalue is the ground-state energy. As e approaches — 1 + , 
the ground-state energy diverges. For e < — 1 there are no 
real eigenvalues. 

There are several ways to obtain the spectrum that 
is displayed in Fig. O. The simplest and most direct 
technique is to integrate the differential equation using 
Runga-Kutta. To do so, we convert the complex differen- 
tial equation fl3.l| ) to a system of coupled, real, second- 
order equations. We find that the convergence is most 
rapid when we integrate along anti-Stokes lines and then 
patch the two solutions together at the origin. This pro- 
cedure, which is described in Ref. gives highly accu- 
rate numerical results. 

To verify the Runge-Kutta approach, we have solved 
the differential equation ( |3.1| ) using an independent and 
alternative procedure. We construct a matrix representa- 
tion of the Hamiltonian in Eq. ([hi]) in harmonic oscillator 
basis functions e~ x / 2 -ff„(a;)7r _1 / 4 j\Jl n n\\ 



M m . n = - 



d.r 



1 



-x 2 /2 



H m (x) 



x e 



dx 2 
- x2/2 H n (x) 



m 



»)] M 2+e } 



(3.2) 



The K-th approximant to the spectrum comes from di- 
agonalizing a truncated version of this matrix M mj „ 
(0 < m, n < K). One drawback of this method is that 



the eigenvalues of M m>n approximate those of the Hamil- 
tonian H in ( fLl| ) only if — 1 < e < 2. Another drawback 
is that the convergence to the exact eigenvalues is slow 
and not monotone because the Hamiltonian H is not Her- 
mitian in a conventional sense. We illustrate the conver- 
gence of this truncation and diagonalization procedure 
for e = -I in Fig. O. 
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FIG. 12. Real eigenvalues of the (K + 1) X (K + 1) 
truncated matrix M m ,n in Eq. (3.2) (K — 0, 1, . . . , 17) for 
e = —1/2. As K increases, the three lowest eigenvalues con- 
verge to the three real energy levels in Fig. hi] at e = —1/2. 
The other real eigenvalues do not stabilize, and instead dis- 
appear in pairs. 

A third method for finding the eigenvalues in Fig. [ll] is 
to use WKB. Complex WKB theory (see Eq. |5~f| ) gives 
an excellent analytical approximation to the spectrum. 

In the next two subsections we examine two aspects of 
the spectrum in Fig. First, we study the asymptotic 
behavior of the ground-state energy as e — > — 1 . Second, 
we examine the phase transition in the vicinity of e = 0. 



A. Behavior of the ground-state energy near e = — 1 

In this subsection we give an analytic derivation of 
the behavior of the lowest real energy level in Fig. [ll] as 
e — > — 1. We show that in this limit the eigenvalue grows 
logarithmically. 



When e = — 1, the differential equation (3.1) reduces 



to 



ip"(x) — ixip(x) = Ei(j(x), 



(3.3) 



which can be solved exactly in terms of Airy functions 
]T^| . The anti-Stokes lines at e = —1 lie at 30° and at 
—210° in the complex- a; plane. We find the solution that 
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vanishes exponentially along each of these rays and then 
rotate back to the real- a; axis to obtain 



lfc,R(aO = Cl,r Ai (rxe^/ 6 + Ee ±2 ™ 



/3 



(3.4) 



We must patch these solutions together at x = accord- 
ing to the patching condition 



-fl^)| 2 
ax 



= 0. 



(3.5) 



But for real E, the Wronskian identity for the Airy func- 
tion Q is 



d_ 

dx 



Ai { xe~™ /6 + Ee- 2in/l 



x=0 



instead of 0. Hence, there is no real eigenvalue. 

Next, we perform an asymptotic analysis for e = —1 + 5 
where 6 is small and positive: 

-ip"(x) - (ix) 1+s *ip(x) = Eip(x), 

1>(x) ~y (x) + 6yi(x) + O(5 2 ) (*->()+). (3.7) 

We assume that E — > oo as 5 — ► 0+ and obtain 

y'o(x) + ixy (x) + Ey (x) = 0, 

y"(x) + ixyi(x) + Eyi(x) = -ix\n(ix)y (x), (3.8) 
and so on. 

To leading order we again obtain the Airy equa- 



tion (3.3) for yo(x). The solution for yo(x) (x > 0) is 
given by i/?r(x) in Eq. (3.4) and we are free to choose 
Cr = 1. We can expand the Airy function in yo(x) for 
large argument in the limit £7 — > oo : 



yo 



-1/4 



(x) = Ai (xe-" /6 + Ee- 2l7! '^ 

~ fxe'^ 6 + Ee~ 2l7r/i 

- (xe~" /6 + Ee- 2i7 '/ i 
3 V 



x exp 



3/2 



(3.9) 



At x = we get 

y (0) = Ai(Ee- 2 ™/ 3 ) 

„ e ^/ e E~ 1 ^e^ E " /2 /(2V^). 



(3.10) 



To next order in e we simplify the differential equation 
for 2/1 (x) in fl3.8|) by substituting 



yi (x) = Q(x)y (x). (3.11) 
Using the differential equation for yo(x) in (|3.8|), we get 
y a {x)Q"(x) + 2y' (x)Q'(x) = -ix\n(ix)y (x). (3.12) 



Multiplying this equation by the integrating factor yo (x) , 
we obtain 



[vl{x)Q' \x)] = -ix\n(ix)y 2 (x), 
which integrates to 



Q'(x) 



Vo( x ) Jx 



dtt ln(it)y 2 (t), 



(3.13) 



(3.14) 



where the upper limit of the integral ensures that Q'(x) 
is bounded for x — > oo. Thus, we obtain 



Q'(o) - 



dxx \a{ix)y Q {x). 



(3.15) 



% 2 (0) ,o 

To determine the asymptotic behavior of the ground- 
state eigenvalue as S — > 0, we insert 

i>{x) ~ y {x) + S yi (x) + 0{S 2 ) 
= y (x) [1 + 5Q(x)] + 0(<5 2 ) 

into the quantization condition: 

x=0 



(3.16) 







dx 
d 

dx 
d_ 
dx 



[\y (x)\ 2 (l + SQ*(x))(l + 5Q(x))} 



x=0 



'M*)l s 



-2% (0)| 2 Re[Q'(0)] 



+25-^ [\y (x)f 



x=0 



Re [Q(0)] 



(3.17) 



We are free to choose Q(0) = 0, and doing so eliminates 
the last term on the right side. The leading-order result 
for the quantization condition in Eq. (|3.6|) then gives 



25|y (0)| 2 Re[g'(0)] 

Z7T 



(3.18) 



Next, we substitute the asymptotic form for yo in 
Eq. (|3.10D and the result for Q'(0) in Eq. fl3.15| ) and ob- 
tain 



2<SRe 



dx ix la(ix) 



Vo(x) 



2/o (0) 



(3.19) 



Because the ratio of t he un perturbed wave functions 
in the integrand in Eq. (3.19) is bounded and vanishes 
exponentially for large x, we know that the integral can 
grow at most as a power of E. Thus, 



CE a e^ E3/2 



(3.20) 



for some power a and constant C and the controlling 
behavior of the ground-state energy as 6 — > is given by 



E 



--\n8 
4 



(3.21) 



where we have neglected terms that vary at most like 
ln(lnJ). Equation (3.21) gives the asymptotic behavior 
of the lowest energy level and is the result that we have 
sought. This asymptotic behavior is verified numerically 
in Table | 
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B. Behavior of energy levels near e = 

In this subsection we examine analytically the phase 
transition that occurs at e = 0. In particular, we study 
high-lying eigenvalues for small negative values of e and 
verify that adjacent pairs of eigenvalues pinch off and 
become complex. 



For small e we approximate H in Eq. (1.1) to first order 



H = p 2 + x 2 + ex 2 ]n(ix) + 0(e 2 ). 



(3.22) 



Using the identity ln(«x) = ha(|x|) + ^iirsgn(x), we then 
have 



H = p 2 



ln (M) + Trsgn(a;) 



+ 0(e 2 ). (3.23) 



The simplest way to continue is to truncate this ap- 
proximate Hamiltonian to a 2 x 2 matrix. We introduce 
a harmonic oscillator basis as follows: The nth eigen- 
value of the harmonic oscillator Hamiltonian p 2 + x 2 is 
E n — 2n + 1 and the corresponding x-space normalized 
eigenstate |n) is 



y/2 n n\ 



(3.24) 



where H n (x) is the nth Hermite polynomial [H (x) 
1 . H^x) = 2x, H 2 {x) = Ax 2 - 2, H 3 (x) = 8x 3 - I2x, 
and so on]. We then have the following diagonal matrix 
elements: 



(n\p 2 + x 2 \n) =271+1, 



(n\x 2 ln(|x|)|n) = a n 
where 7 is Euler's constant and 



l +ln2 ) ( n+1 2 



[n+l/2] 



l + [n/2] + (n + l/2) £ 



1 



2k- 1 



(3.25) 



(3.26) 



(3.27) 



We also have the off-diagonal matrix element 

TABLE I. Comparison of the exact ground- state energy E 
near e = —1 and the asymptotic results in Eq. (3.21). The ex- 
plicit dependence of E on e = —1+5 is roughly E oc (— In S) 2 ^ 3 
as 5 -> 0+. 



E c 



Eq. (3.21) 



0.1 

0.01 

0.001 

0.0001 

0.00001 

0.000001 

0.0000001 



1.6837 
2.6797 
3.4947 
4.1753 
4.7798 
5.3383 
5.8943 



2.0955 
2.9624 
3.6723 
4.3013 
4.8776 
5.4158 
5.9244 



(2n — l|-i7ra; 2 sgn(a:)|2n) 



~i(8n+l) 



r 2 (n+l/2) 



n\(n- 1)! 



1/2 



(3.28) 



In the (2n — 1) — (2n) subspace, the matrix H ~E then 
reduces to the following 2x2 matrix: 



A-E iB 
iB C-E 



where for large n and small e we have 



(3.29) 



A ~ An - 1 + e(n - 1/2) ln(2n), 
g 

B ~ —en, 
3 ' 

C ~ 4n+ 1 + enln(2n). 



(3.30) 



The determinant of the matrix in Eq. ( 3.2S ) gives the 
following roots for E: 



1 



E= - [A + C±y/{A-C) 2 -4B 2 



(3.31) 



We observe that the roots E are degenerate when the 



discriminant (the square-root) in Eq. (3.31) vanishes. 
This happens when the condition 



3 

8n 



(3.32) 



is met. Hence, the sequence of points in Fig. [ll] where 
the eigenvalues pi nch o ff approaches e = as n — > 00. 
For example, Eq. ( 3.3 2| ) predicts (using n = 4) that E7 
and Eg become degenerate and move off into the complex 
plane at e w — 0.1. In Fig. [l^ we compare our prediction 
for the behavior of E in Eq. ( 3.31 ) with a blow-up of 
of a small portion of Fig. [ll]. We find that while our 
prediction is qualitatively good, the numerical accuracy 
is not particularly good. The lack of accuracy is not 
associated with truncating the expansion in powers of e 
but rather with truncating the Hamiltonian H to a 2 x 2 
matrix. Our numerical studies indicate that as the size of 
the matrix truncation increases, we obtain more accurate 
approximations to the behavior of the energy levels E in 

Fig. 

IV. MORE GENERAL CLASSES OF THEORIES 

In this section we generalize the results of Sees. [0] and 



I [J to a much wider class of theories. In particular, we 
consider a complex deformation of the x 2K anharmonic 
oscillator, where K = 1, 2, 3, ... [see Eq. (|lJ)]. The 
Schrodinger equation for this oscillator has the form 



tl>"(x) + [x 2K (ix) e - E]i/)(x) = 0. 



(4.1) 
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FIG. 13. A comparison of the prediction in Eq. (3.31) 
and a magnification of Fig. [Tl| Our prediction for the point 
at which E7 and Es become degenerate is not very accurate 
numerically but is qualitatively quite good. 



To determine the energy levels E as functions of the 
deformation parameter e, we must impose appropriate 
boundary conditions on Eq. fl4.l| ). We require that the 
wave function vanish as |x| — * 00 inside of two wedges 
symmetrically placed about the imaginary- a; axis. The 
right wedge is centered about the angle bright j where 



C7T 



'right 



4K + 2e + 4' 



(4.2) 



and the left wedge is centered about the angle #i e ft, where 

C7T 



0\eit = — 7T ' 



4K + 2e + 4 
The opening angle of each of these wedges is 
2tt 



2K + e + 2 



(4.3) 



(4.4) 



This pair of wedges is PT (left-right) symmetric. 

The orientation of these wedges is determined by an- 
alytically continuing the differential equation eigenvalue 



problem (4.1) and associated boundary conditions in the 
variable e using the techniques explained in Ref. ]lq ] . The 
rotation of the boundary conditions is obtained from the 
asymptotic behavior of the solution ip{x) for large |x|: 



ip{x) ~ exp I ± 



je/2 a ,_R'+l+e/2 

K +l + e/2 



(4.5) 



(In this formula we give the controlling factor of the 
asymptotic behavior of the wave function; we neglect al- 
gebraic contributions.) Note that at the center of the 
wedges the behavior of the wave function is most strongly 



exponential; the centerline of each wedge is an anti- 
Stokes line. At the edges of the wedges the asymptotic 
behavior is oscillatory. The lines marking the edges of 
the wedges are Stokes lines. 

For all positive integer values of K the results are qual- 
itatively similar. At e = the two wedges are centered 
about the positive and negative real axes. As e increases 
from the wedges rotate downward and become thinner. 
In the region e > the eigenvalues are all real and posi- 
tive and they rise with increasing e. As e — > 00, the two 
wedges become infinitely thin and lie along the negative 
imaginary axis. There is no eigenvalue problem in this 
limit bec ause the solution contour for the Schrodinger 
equation (4.1) can be pushed off to infinity. Indeed, we 
find that in this limit the eigenvalues all become infinite. 

When e is negative, the wedges rotate upward and be- 
come thicker. The eigenvalues gradually pair off and 
become complex starting with the highest eigenvalues. 
Thus, VT symmetry is spontaneously broken for e < 0. 
Eventually, as e approaches —K 1 only the lowest eigen- 
value remains real. At e = —K the two wedges join at the 
positive imaginary axis. Thus, again there is no eigen- 
value problem and there are no eigenvalues at all. In the 
limit e — > — K the one remaining real eigenvalue diverges 
logarithmically. 

The spectrum for the case of arbitrary positive integer 
K is quite similar to that for K = 1. However, in general, 
when K > 1, a novel feature emerges: A new transition 
appears for all negative integer values of e between and 
—K. At these isolated points the spectrum is entirely 
real. Just above each of these negative-integer values of 
e the energy levels reemerge in pairs from the complex 
plane and just below these special values of e the energy 
levels once again pinch off and become complex. 



A. Quantum x 4 (ix) e theory 

The spectrum for the case K = 2 is displayed in 
Fig. |l4]. This figure resembles Fig. [ll] for the case K = 1. 
However, at e = — 1 there is a new transition. This tran- 
sition is examined in detail in Fig. 

An important feature of the spectrum in Fig. 
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the disappearance of the eigenvalues and divergence of 
the lowest eigenval ue as e decreases to —2. Following 
the approach of Sec. Ill A, we now derive the asymptotic 
behavior of the ground-state energy as e ^ — 2 + . To 
do so we let e = — 2 + 5 and obtain from Eq. (4.1) the 
Schrodinger equation 



ip"(x) — x 2 (ix) s i/j(x) = Etj}(x). 



(4.6) 



We study this differential equation for small positive 5. 
When 6 = this differential equation (4.6) reduces to 



i>"(x) - x 2 ip(x) = Eip(x). 



(4.7) 



13 



50 



40 



dx 

L 

a 



30 



20 



10 




FIG. 14. Energy levels of the Hamiltonian 

as a function of the parameter e. This figure 
is similar to Fig. [H], but now there are four regions: When 
e > 0, the spectrum is real and positive and it rises mono- 
tonically with increasing e. The lower bound e = of this 
"PT-symmetric region corresponds to the pure quartic anhar- 
monic oscillator, whose Hamiltonian is given by H = p 2 + 
When — 1 < e < 0, VT symmetry is spontaneously broken. 
There are a finite number of real positive eigenvalues and an 
infinite number of complex conjugate pairs of eigenvalues; as 
a function of e the eigenvalues pinch off in pairs and move off 
into the complex plane. By the time e = — 1 only eight real 
eigenvalues remain; these eigenvalues are continuous at e = 1. 
Just as e approaches —1 the entire spectrum reemerges from 
the complex plane and becomes real. (Note that at e = — 1 
the entire spectrum agrees with the entire spectrum in Fig. |ll| 
at e = 1.) This reemergence is difficult to see in this figure 
but is much clearer in Fig. [l^ in which the vicinity of e = — 1 
is blown up. Just below e = — 1, the eigenvalues once again 
begin to pinch off and disappear in pairs into the complex 
plane. However, this pairing is different from the pairing in 
the region — 1 < e < 0. Above e = — 1 the lower member of a 
pinching pair is even and the upper member is odd (that is, 
E& and Eg combine, E\o and En combine, and so on); below 
e = — 1 this pattern reverses (that is, E-j combines with E$, 
Eg combines with E\o, and so on). As e decreases from —1 to 
—2, the number of real eigenvalues continues to decrease until 
the only real eigenvalue is the ground-state energy. Then, as 
e approaches — 2 + , the ground-state energy diverges logarith- 
mically. For e < —2 there are no real eigenvalues. 

The anti-Stokes lines for this equation lie at 45° and at 
—225°. Thus, we rotate the integration contour from the 
real axis to the anti-Stokes lines and substitute 



V2' 



V2 



e * 



(Re x < 0) 
(Re x > 0) 



(4.8) 



for x in the left-half and in the right-half complex plane, 
respectively. Note as s and r increase, x moves towards 
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FIG. 15. A magnification of Fig. [ll| in the vicinity of the 
transition at e = —1. Just above e = — 1 the entire spectrum 
reemerges from the complex plane, and just below e = — 1 it 
continues to disappear into the complex plane. The spectrum 
is entirely real at e — — 1. 

complex infinity in both the left- and right-half plane. 

The wave function in the left-half plane, V'l(s), and 
the wave function in the right-half plane, tp^r), satisfy 
the differential equations 



dr 2 



W>l(s), 

Mr) = (-v-l)Mr), (4.9) 



where we have set v = — iE — i. For each of these equa- 
tions the solution that vanishes at infinity is a parabolic 
cylinder function pq] : 

lM«) - C L D v (s) = ClD, (xV2e^ , 

Mr) = C R U_ v _i(r) = C R D_ V _ X (xyfar*) , (4.10) 

where Cl and Cr are arbitrary constants. 

We impose the quantization condition by patching 
these solutions together at x = on the real- a; axis ac- 
cording to the patching conditions 



x=0 



dx 



x=0 

-^r(z) 



x=0 



(4.11) 



d 

x =o dx 

To eliminate the constants Cl and Cr we take the ratio 
of these two equations and simplify the result by cross 
multiplying: 



dx dx 



0. 



(4.12) 



x=0 
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We now show that this condition cannot be satisfied 



by the 6 = wave function in Eq. (4.10). For this case, 
the quantization condition (4.12) states that 



D v (s)4-D-v-i(ia) - D- v -x (is)-^-D u (s) (4.13) 
as as 

vanishes at s = 0. (We have simplified the argument by 
setting s = x \/2e 5iTr/4 .) But Eq. Q4.13| ) for any value of s 
is just the Wronskian for parabolic cylinder functions ]l6| ] 
and this Wronskian equals — ze~'" y7r / 2 . This is a nonzero 
result. Thus, when 8 = 0, there cannot be any eigenvalue 
E, real or complex, and the spectrum is empty. 



The quantization condition (4.12) can be satisfied 
when 6 > 0. We investigate this region for the case 
when 6 is small and positive by performing an asymp- 
totic analysis. We assume that E — > oo as 6 — » 0+, but 
slower than any power of 8, and that the wave function 
ip(x) has a formal power series expansion in 8: 

i>{x)~y (x)+8 yi (x) + O(8 2 ) (8^0+). (4.14) 



Next, we expand the Schrodinger equation (4.6) in pow- 
ers of 8: 

y'J(x)+x 2 y (x)+Ey (x) = 0, 

y'({x) + x 2 yi{x) + E yi (x) = -x 2 ln(ix)y (x), (4.15) 
and so on. 



Of course, to zeroth order in 8 we obtain Eq. (4.7) for 
yo(x). Thus, in the left- and right-half complex x-planc 
we get 



ykix) = C L D„ (xV2e^) , 
y*(x) = C K D_ V _ X (aVar*) • (4.16) 
To first order in 8, we simplify the differential equation 



for yi{x) in (4.15) by substituting 
yi[x) = Q(x)y (x). 



(4.17) 



Using the differential equation for yo(x) in (4.15), we get 

y (x)Q"(x) + 2y' (x)Q'(x) = -x 2 ln(ix)y (x). (4.18) 

Multiplying this equation by the integrating factor yo(x), 
we obtain 



[y 2 (x)Q'(x)]' = -x 2 \n(ix)y 2 (x). 
The integral of this equation gives 



(4.19) 



t/4 



dtt 2 ln(it) 



dtt 2 In (it) 



yb( x ) 
y R (i) ' 



(4.20) 



where the limit of the integral at infinity ensures that 
Q'(x) is bounded for \x\ — > oo. 

To determine the asymptotic behavior of the ground- 
state eigenvalue as 8 — ► + , we insert 



= y^ R (x) [1 + SQ h > R (x)] 



L,R, 



(4.21) 



into the quantization condition ( 4.12j ) 




yo( x )-^yo( x )-yo(x)-^y^(x) 



x=0 



x=0 



x [1 + *(Qr(0) + Ql(0))] 
+^ R (0)2/ L (0) [Ql(0) - Q' R (0)] 



(4.22) 



We are free to choose Qr(0) + Ql(0) = to simplify this 
result. 

Substituting the Wronskian for the parabolic cylinder 
function and the result for yo(0) in Eq. ( 4.16| ), we obtain 



V2e-^ ~ 8D v (p)D- v -i(0) [Q' R {0) - Q' L (0)} . 
We can simplify this result using the identity 



D v {0) 



and v = — \E — ~ to obtain 



£>„(0).D_„-i(0) 



r (i + 1) cos 



H) 



e 4 



r(i + |)V2 



(4.23) 



(4.24) 



(4.25) 



where we have used the reflection formula, T(z)T(l — z) = 
it/ sin(7rz), and the asymptotic behavior r(x+l/2)/r(x+ 
1) ~ x" 1 ^ 2 for large x. Thus, Eq. ( 4.23j ) reduces to 



'2E 



Qk(o) - QL(o). 



(4.26) 



We can further show that 



Qk(o)-OL(o) 

dtt 2 In(it) 



LU-iC^te"^) 



D. 



v-l 



(0) 



dtt 2 ln(it) 



(i 



D v {y/2te***/*) 



°°t 2 dt^, ( S 3^ 

— 5— e 4 In — — e 4 

o 2§ \V2 



D v {0) 



D. 



t(*) 



°°t 2 dt j,., 

— s— e 4 in I — =e 
o 2§ VV2 



(JL t 



D v {t) 



D u {0) 



(4.27) 
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We observe that the previous expression is real because 
v* = implies that D v (t)* = D- v -i(t) and t hus t he 

two integrals are complex conjugates. Thus, Eq. ( 4.27 ) is 
real, and E is a real function of 6. Furthermore, because 
the ratio D v {t)/D V {Q) appears in both integrals, the ex- 
pression can at most vary as a power of E. Hence, the 
contribution of Q' R (0)-Q' L (0) to the balance in Eq. fl4.26j ) 
is subdominant and can be neglected. Our final result for 
the small-5 behavior of the lowest eigenvalue is that 



E ln5 + 0[ln(ln(5)] (5->0+). 



(4.28) 



In Fig. we show that Eq. ( 4.28| ) compares well with 
the numerical data for the lowest eigenvalue in the limit 
as 5 — > 0. 
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FIG. 16. A comparison of the lowest eigenvalue of the 
Hamiltonian H = p 2 + x 4 (ixY (solid circles) with the asymp- 
totic prediction in (4.28) (solid line) near e = —2. The 
solid line includes a one parameter fit of terms that grow like 
ln(ln<5) as S -» 0+. 



B. Classical x (ixY theory 

It is instructive to compare the quantum mechanical 
and classical mechanical theories for the case K = 2. Our 
objective in doing so is to understand more deeply the 
breaking of VT symmetry that occurs at e = 0. For the 
case K = 1 we found that VT symmetry is broken at the 
classical level in a rather obvious way: Left-right sym- 
metric classical trajectories become spirals as e becomes 
negative (see Fig. 9). However, we find that when K = 2 
spirals do not occur until e < — 2. The classical manifes- 
tation of VT symmetry breaking for — 2 < e < and the 
transition that occurs at e = is actually quite subtle. 

For purposes of comparison we begin by examining the 
classical trajectories for the positive value e = 0.7. In 



Fig. 17 we plot three classical trajectories in the complex- 
x plane. The first is an arc that joins the classical turning 
points in the lower-half plane. The other two are closed 
orbits that surround this arc. The smaller closed orbit re- 
mains on the principal sheet and has a period (T « 4.9), 
which is equal to that of the arc. The more complicated 
trajectory is left-right symmetric but extends to three 
sheets of the Riemann surface. The period (T 26.1) of 
this third orbit is significantly different from and larger 
than the period of the other two. 



FIG. 17. Three classical trajectories in the complex- a; plane 
for a particle described by the Hamiltonian H = p 2 + x 4 (ixY 
with e = 0.7. The solid line represents oscillatory motion 
between the classical turning points. The long-dashed line is 
a nearby trajectory that encloses and has the same period as 
the solid-line trajectory. The short dashed line has a different 
topology (it enters three sheets of the Riemann surface) from 
the long-dashed line, even though these trajectories are very 
near one another in the vicinity of the turning points. The 
period of this motion is much longer than that of the solid 
and long-dashed trajectories. 

Next, we consider the negative value e = — 0.7. In 
Fig. |lj we plot two classical trajectories for this value. 
The first (solid line) is an arc joining the classical turn- 
ing points in the upper-half plane. This arc extends to 
three sheets of the Riemann surface. The other trajec- 
tory (dashed line) is a closed orbit that surrounds this 
arc. Both have the period T « 22.3. This figure illus- 
trates the first of two important changes that occur as e 
goes below zero. The trajectory that joins the two turn- 
ing points no longer lies on the principal sheet of the 
Riemann surface; it exhibits a multisheeted structure. 

Figure |l^ illustrates the second important change that 
occurs as e goes below zero. On this figure we again plot 
two classical trajectories for the negative value e = —0.7. 
The first (solid line) is the arc joining the classical turning 
points in the upper-half plane. This arc is also shown on 
Fig. EsL The second trajectory (dashed line) is a closed 
orbit that passes near the turning points. The two trajec- 
tories do not cross; the apparent points of intersection are 
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on different sheets of the Riemann surface. The period of 
the dashed trajectory is T w 13.7, which is considerably 
smaller than that of the solid line. Indeed, on the basis 
of extensive numerical studies, it appears that all tra- 
jectories for — 2 < e < 0, while they are VT (left-right) 
symmetric, have periods that are less than or equal to 
that of the solid line. When e > 0, the periods of tra- 
jectories increase as the trajectories move away from the 
oscillatory trajectory connecting the turning points. 

8 i 1 1 1 1 1 1 1 1 
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FIG. 18. Two classical trajectories in the complex- a; plane 
for a particle described by the Hamiltonian H = p' 2 + x 4 (ixY 
with e = —0.7. The solid line represents oscillatory motion 
between the classical turning points. This trajectory enters 
three sheets of the Riemann surface. The dashed line is a 
nearby trajectory that encloses and has the same period as 
the solid-line trajectory. 
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FIG. 19. Two classical trajectories in the complex- a; plane 
for a particle described by the Hamiltonian H = p 2 + x 4 (ix) e 
with e = —0.7. The solid line represents oscillatory motion 
between the classical turning points and is the same as that in 
Fig. 18. The dashed line is a nearby trajectory whose period 
is smaller than the period of the solid-line trajectory. 

We speculate that for negative values of e the appear- 
ance of complex eigenvalues in the quantum theory (see 



Fig. 14) is associated with an instability. The path in- 
tegral for a quantum theory is ordinarily dominated by 
paths in the vicinity of the classical trajectory connect- 
ing the turning points. However, when e is negative, we 
believe that these trajectories no longer dominate the 
path integral because there are more remote trajectories 
whose classical periods are smaller. Thus, the action is 
no longer dominated by a stationary point in the form of 
a classical path having VT symmetry. Hence, the spec- 
trum can contain complex eigenvalues. 

The appearance of a purely real spectrum for the spe- 
cial value e = — 1 is consistent with this conjecture. For 
integer values of e > —2 we hnd that all classical trajec- 
tories lie on the principal sheet of the Riemann surface 
and have the same period. 

C. Quantum x 6 (ix) e theory 

The spectrum for the case K = 3 is displayed in 
Fig. ^0|. This figure resembles Fig. |lj for the case K = 2. 
However, now there are transitions at both e = — 1 and 
e = -2. 

V. COMPLEX DEFORMATIONS OF 
NONANALYTIC POTENTIALS 

In our discussion so far we have considered complex 
deformations of the potentials x 2K . These potentials 
are analytic functions of x. In this section we consider 
complex deformations of the nonanalytic potentials \x\ p , 
where P is real. We will see that the eigenvalues of the 
potential |a;| p (ia;) e are real only when e = (and some- 
times at other isolated values of e). Thus, it appears that 
if one attempts to construct a complex deformation of a 
nonanalytic potential, one destroys a crucial property of 
the theory; namely, that the spectrum be real. 

We begin our discussion by examining the spectrum 
of the \x\ p potential. The eigenvalues of this potential 
are displayed as a function of P in Fig. |2l]. It is inter- 
esting that the spectrum of this potential Jl7j is quite 
similar to that of the x 2 (ix) e potential for positive e (see 
Fig. |llj) . The difference between the spectra of these two 
potentials becomes apparent when e is large: As e — > oo, 
the spectrum of |x| 2+c approaches that of the square- well 
potential [E n = (n + l) 2 7r 2 /4], while the energies of the 
x 2 (ix) e potential diverge. 

WKB theory gives an excellent approximation to the 
spectrum of both potentials and thus provides an inter- 
esting comparison. For the x 2 (ix) e potential, when e > 0, 
the novelty of the WKB calculation is that it must be per- 
formed in the complex plane. The turning points x± are 
those roots of E — x 2 (ix) e = that analytically continue 
off the real axis as e moves away from zero (the harmonic 
oscillator): 
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= E— e' 



4 + 3c 



(5.1) 



These turning points lie in the lower-half (upper-half) x 
plane when e > (e < 0). 



r(f±ff) y/Z(n + l/2) 
_ sin(^)r(|±|) 



(n 



oo 



(5.4) 





FIG. 20. Energy levels 

of the Hamiltonian H — p 2 + x 6 (ix) £ as a function of the 
parameter e. This figure is similar to Fig. |l4|, but now there 
are five regions: When e > 0, the spectrum is real and pos- 
itive and it rises monotonically with increasing e. The lower 
bound e = of this "PT-symmetric region corresponds to 
the pure sextic anharmonic oscillator, whose Hamiltonian is 
given by H — p 2 +x 6 . The other four regions are —1 < e < 0, 
-2 < e < -1, -3 < e < -2, and e < -3. The PT symmetry 
is spontaneously broken when e is negative, and the num- 
ber of real eigenvalues decreases as e becomes more negative. 
However, at the boundaries e = — 1, — 2 there is a complete 
real positive spectrum. When e = — 1, the eigenspectrum is 
identical to the eigenspectrum in Fig. h4at e = 1. For e < —3 
there are no real eigenvalues. 

The leading-order WKB phase-integral quantization 
condition is 



2n+ 1 



dx y/E — x 2 {ixY 



(5.2) 



It is crucial that this integral follow a path along which 
the integral is real. When e > 0, this path lies entirely in 
the lower-half x plane and when e = the path lies on 
the real axis. But, when e < the path is in the upper- 
half x plane; it crosses the cut on the positive-imaginary 
axis and thus is not a continuous path joining the turning 
points. Hence, WKB fails when e < 0. 

When e > 0, we deform the phase-integral contour so 
that it follows the rays from x~ to and from to x+: 



2n + l 



tt = 2 sin (^j— ) E ^ J o <Wl-s 2+e . (5.3) 



We then solve for E n : 



FIG. 21. Energy levels of the Hamiltonian H — p 2 + \x\ p 
as a function of the parameter P. This figure is similar to 
Fig. |ll], but the eigenvalues do not pinch off and go into the 
complex plane because the Hamiltonian is Hermitian (spec- 
trum becomes dense at P — 0). 

We can perform a higher-order WKB calculation by 
replacing the phase integral by a closed contour that en- 
circles the path connecting the two turning points (see 
Ref. 111). With Q(x) = x 2 (ix) e -E, the next to leading- 
order WKB quantization condition is 



2n+ 1 



1 



2i 



7T = — f dx y/Q(x) + — 



2 1 



, 48Q(a;)5 



(5.5) 



where the contour C encircles the turning points x+ and 
X- in a counterclockwise direction. As above, we deform 
the contour to lie above and below the rays that connect 
the turning points with x = 0, we obtain for the second 



contour integral in Eq. (5.5): 



lldx^l 
2iJc 48Q(x)5 

_jt±^2 + 3e-f 

= —E 4 + 2 « 
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4+ g e 
-E — sin , 
24 V2 



■ sin 



ds — 
" (1 



r ( 3+ 2e 

< ) v ( 5 + 3e 
1 U+ 2e 



s 2+e) 2 

(5.6) 



We add the contribution of this integral to the leading- 
order result in Eq. (5.3) and solve the resulting expression 
for the energy E, assuming that n is large, and obtain 
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En 



r ( 8+3e 
1 U+ 2e 



) 0r(n + 1/2) 



sin(^)r(f±f) 

(2 + e)(l + e)sin(^) 
6^ (n+i) 2 (4 + e) 2 



(n->oo). (5.7) 



This is the next-to-leading-order WKB result for the en- 
ergy. 

The next correction to the WKB result for the |a;| 2+e 
potential is 



E n 



1 - 2r 
4 + e 



r (§rl) V^(" + i/2) 

(2 + e)(l + e)cot(^) 
3tt (n+i) 2 (4 + e) 2 



(n — * oo), (5. 



Note that the leading-order WKB quantization condi- 
tion (accurate for e > —2) for the |x| 2+e potential is like 



Eq. (5.4) except that sin (^2TeJ * s absent. 
10 
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FIG. 22. Energy levels 

of the Hamiltonian H = p 2 + \x\(ix) e as a function of the 
parameter e. The spectrum is entirely real only when e = 0. 



into the complex- a; plane; we cannot unambiguously de- 
fine the rotation of wedges in which the boundary condi- 
tions apply. However, for sufficiently small e we can allow 
x to remain real and we can impose the boundary con- 
dition that ip(x) — ► for x — * ±oo. Specifically, we have 
the condition for all P that if |e| < 2, then this boundary 
condition on the real- a: axis may be consistently imposed 
to define the eigenspectrum. 

We now consider two cases: P = 1 (Fig. 22) and P = 3 
(Fig. |23| ). Fig ure |22| is quite similar to Fig. 11 and Fig. |23| 
resembles Fig. |14|. The key properties of these figures are 
that (1) the lowest energy level diverges at e = — P/2, 
and that (2) the energy levels pinch off and go into the 
complex plane on both sides of e = 0. Thus, the spectrum 
is entirely real only when e = 0. 




FIG. 23. Energy levels as a function of the parameter e 
for the Hamiltonian H = p 2 + |a;| 3 (ia;) e . The spectrum is real 
when e = and e = —0.5. 
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Now we examine what happens when we attempt to 
deform the \x\ p potential using a complex deformation. 
That is, we consider an |a;| p (ia;) e potential. Of course, 
since \x\ is not an analytic function, we cannot define 
an analytic continuation of the Schrodinger eigenvalue 
problem 

d 2 

~ T^( x ) + \x\ p (txY^{x) = Ei/j(x) (5.9) 
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